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ABSTRACT 

Linear stability of a current sheet that is subject to an impulsive acceleration due to a shock passage is studied 
with the effect of guide magnetic field. We find that the current sheet embedded in relativistically magnetized 
plasma always shows a Richtmyer-Meshkov type instability, while it depends on the density structure in the 
Newtonian limit. The growth of the instability is expected to generate turbulence around the current sheet that 
can induce so-called turbulent reconnection whose rate is essentially free from plasma resistivity. Thus, the 
instability can be applied as a triggering mechanism of rapid magnetic energy release in variety of high-energy 
astrophysical phenomena such as pulsar wind nebulae, gamma-ray bursts, and active galactic nuclei, where the 
shock wave is supposed to play a crucial role. 

Subject headings: instabilities — magnetic fields — shock waves — relativistic processes — turbulence 



1. INTRODUCTION 

Release of magnetic energy is believed to be essentially 
important in high-energy astrophysical phenomena such as 
pulsar wind nebulae, gamma-ray bursts, and active galactic 
nuclei (e.g., Kennel & Coroniti 1984, Giannios et al. 2009, 
Lyubarsky 2010). However, the detailed physics of the dissi- 
pation mechanism of the magnetic field is still unclear. Theo- 
retical studies have been suggesting that instability of current 
sheet or turbulent environment enhances the rate of magnetic 
reconnection substantially independent of plasma resistivity 
(Lazarian & Vishniac 1999, Kowal et al. 2009). Takamoto 
et al. (2012) also showed that turbulent stretching of current 
sheet leads to the dissipation of the magnetic field in turbulent 
eddy quickly in a few eddy turnover time regardless of plasma 
resistivity. 

Recently, the current sheet embedded in relativistically 
magnetized plasma under the effect of secular acceleration 
was found to be unstable (Lyubarsky 2010). Here the rela- 
tivistically magnetized plasma means that the magnetic en- 
ergy density is larger than rest mass energy density. De- 
spite the instability found by Lyubarsky is a Rayleigh-Taylor 
(or Kruskal-Schwarzchild) type instability, it can be unstable 
even when the system is isochoric, since the magnetic and the 
thermal energies play the role of inertia in the relativity. The 
Richtmyer-Meshkov instability, which is induced by an im- 
pulsive acceleration due to shock passage (Richtmyer 1960, 
Nishihara et al. 2010), is a counterpart of the Rayleigh-Taylor 
instability, suggesting that the shock-current sheet interaction 
in the relativistically magnetized plasma can be unstable as a 
counterpart of Lyubarsky's instability. Since the shock wave 
is an essential ingredient in the high-energy astrophysical phe- 
nomena, the shock-current sheet interaction is quite ubiqui- 
tous. 

For these reasons, we study linear stability of the relativistic 
current sheet that is subject to an impulsive acceleration due 
to the shock passage. The organization of this paper is as 
follows. In §2, we provide unperturbed zeroth-order state of 
the current sheet. Then, in §3, the linear stability analysis 
is performed and the master equation that governs instability 
growth is derived. In §4, basic properties of the instability is 
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studied based on the solution of the master equation. Finally, 
in §5, we discuss implications of the instability. 

2. ZEROTH-ORDER STATE 

In this paper, we use units of the speed of light c = 1 and 
the magnetic permeability /i = 1 . We consider a static initial 
current sheet as follows: Magnetized media that have oppo- 
sitely oriented z-component magnetic field are separated by a 
current sheet. In the left and right sides of the current sheet 
(henceforth region 1 and region 3, respectively), the magnetic 
field is in y-z plane Z?i(3) = (0,-8-^1(3), fi Zl i(3)), where B-\ =—B z ^. 
The y-component magnetic field, so called the guide field, 
is constant across the current sheet (B y \ = B y 2 = B y ^ = B y ). 
Inside the current sheet, where the z-component magnetic 
field is dissipated, the total pressure is balanced with ex- 
ternal plasma as /?2 + fi 2 2 /2 = /?i ( 3) + Z? 2 1(3) /2 + .B 2 1(3) /2. If 
the plasma in region 1 (and 3) are relativistically magnetized 
and cold (B Zi \ p\,p\), the ratio of the inertia (or total en- 
thalpy) of the media 1 and 2 is w 2 /wi = {p2 + 72^2/(72- 1) + 
^ 2 2 }/{Pi + 7iPi/(7i-l)+4i+^.i}-(2^. 1 +3^ 2 1 )/(fi z 2 1 + 
B 2 j), where the adiabatic index 72 =4/3 is used in the last 
expression for the relativistically hot plasma of the region 2 
(pi <C P2 ~ 0[B\~\). In particular, when there is no guide 
field, the ratio becomes W2/W1 ~ 2. This indicates that the 
current sheet is heavier than the external plasma, and this is 
the essence of the instability shown below. 

Let us consider the propagation of a fast relativistic- 
magnetohydrodynamic (RMHD) shock wave in the region 1 
toward +x direction. In the relativistically magnetized cold 
plasma (B\ > p\ > p\), the shock velocity is close to the 
speed of light (v s h ~ 1), because the fast characteristic speed 
is close to the light velocity, while the fast shock causes only 
a small jump of fluid velocity (also density and magnetic field 
strength as well). According to the jump condition of the 
RMHD shock for such plasma (Kennel & Coroniti 1984), the 
post shock plasma velocity in the upstream rest-frame is given 
by 

V ^sh-2(1-Vsh)fi 2 /Pi 
vl h + 2{\-v^B\/ Px ' 

By solving Eq. (fl}, we see that the post shock velocity 
can be subsonic (V < c s = 1/^/3), provided B\j ' p\ > v 2 h (l - 
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FIG. 1. — Schematic of the system under consideration. Shaded areas are 
the regions for the perturbations that subject to the acceleration compared to 
the zeroth order state. 

CsV s h)/{2(l +c s )(l - v s h)}. Thus, for instance, even in the 
cases of v s h = 0.9, 0.99, and 0.999, the post shock velocity 
is subsonic for the medium of B\/p\ > 1.23, 13.3, and 134, 
respectively. The postshock velocity V basically gives the de- 
gree of impulsive acceleration for the current sheet. We dis- 
cuss about this point in more detail in the next section. 

3. LINEAR STABILITY ANALYSIS 

If the initial current sheet is uniform, the shocked media 
are accelerated and move homogeneously after the shock pas- 
sage. However, when the current sheet is perturbed as illus- 
trated in Fig.Q] the shock acceleration near the current sheet 
acts inhomogeneously depending on the y-coordinate due to 
the different inertia of the media. In the following, we ex- 
amine the dynamics of perturbations in the rest-frame of the 
homogeneous shocked current sheet, indicating that the ze- 
roth order variables are taken as their postshock values when 
the current sheet is uniform. We denote the left (right) side of 
the interface by the interface A (B), its amplitude by £a (£b), 
and the mean position by x = (x = A). 

We consider perturbation due to an impulsive shock passage 
that is imposed as a force term on the basic equations given 
below following the formulation by Wheatley et al. (2005), 
and the linearization is performed in order to find the response 
of such a perturbation. Since the perturbed flow can be treated 
as subsonic one owing to eq. (Q]l, we can assume the Lorentz 
factor of the perturbed flows to be unity and employ the in- 
compressive approximation. Then, if we write the perturbed 
variables by using the small letters and the unperturbed zeroth 
order variables by the capital letter, the basic RMHD equa- 



tions for the perturbations can be written as 

d xVx + dyVy = 0, (2) 

(W+Bl+Bl)d t v x =-d x p-B y d x b y +B y dyb x +F{t,x,y), (3) 

{W + Bl)d,v y = -d y p, (4) 

d t b X =BydyV X , (5) 

d t by=-Byd X V X , (6) 

d X b X + dyby = 0, (7) 



where W = p+jP/Qy— 1) is the specific enthalpy and F(t,x,y) 
is the force term that represents the acceleration of the per- 
turbed flow. In the above expression, we have omitted the 
z-component velocity and magnetic field perturbations owing 
to the fact that they are decoupled from the other variables 



and only describe the propagation of the Alfven wave in y- 
direction. 

Because the shock crossing time of the current sheet can 
be much smaller than the timescale of perturbations, the force 
term can be modeled as 

F(t,x,y) = (W 2 + Bl 2 + Bl 2 -W l -Bl l -Bl l ) 

x{H[x]-H[x-Uy,t)]}V5(t) 

+(W 3 + If,,-, + B 2 z3 -W 2 - B\ 2 - «\) 

x {H[x+A]-H[x+A-^ 3 (y,t)]}V6(t), (8) 

where H(t ) is the Heaviside function and 5{t ) is the Dirac delta 
function (see also Wheatley et al. 2005). The regions for the 
perturbations that subject to the acceleration compared to the 
zeroth order state are illustrated as shaded regions in Fig.[T] In 
this formulation, we have to careful about following points: In 
Wheatley et al. (2005) that studied the Richtmyer-Meshkov 
instability (instability of a single contact surface), the speed 
of contact surface after the shock passage is chosen as V in 
eq. dHJ, which is slightly different from the postshock ve- 
locity. Thus, if we consider the case where there is only the 
interface A, the speed of the contact (interface A) between the 
shocked regions 1 and 2 must be used as V in eq. (0, which 
is obtained by solving an appropriate Riemann problem (see, 
e.g., Giacomazzo & Rezzolla 2006). This suggests that for the 
stability of the interface A in the short-wavelength limit (see 
eq. |29l below), the substitution of the speed of the shocked 
contact as V in eq. (O rather than the postshock velocity of 
eq. ([TJ would give more accurate growth rate. For the fi- 
nite thickness current sheet, the speed of the interfaces A and 
B after the shock passage approaches the postshock velocity 
V of eq. ([T} asymptotically with time. Thus, since the shock 
crossing time of the current sheet is much smaller than growth 
timescale of the instability, the post shock velocity of eq. (fl3 
is more appropriate as the degree of the impulsive accelera- 
tion V in eq. (HJ. Another point we have to keep in mind is 
the neglect of the effect of a reflection shock (or a rarefac- 
tion wave) that is generated when the incident shock hits the 
interface B, which causes complex dynamics in the current 
sheet. However, direct numerical simulations of the (Newto- 
nian and unmagnetized) Richtmyer-Meshkov instability of a 
finite thickness fluid layer showed that the impulsive model 
can indeed work well, although it cause error about a factor 2 
or less in growth rate (Mikaelian 1996). 

Assuming the perturbations of the form q(x,y,t) = 
q(x,t) exp(iky) and taking the temporal Laplace transform 
(q[x,s] = L q[x,t] exp[-sf] dt) of ©-([Til outside of the 
forced region give 

d x v x +ikv y = 0, (9) 
s ( W + B] + B z z ) v x = -d x p- B y d x b y + ikB y b x , (10) 
s(W+B 2 z )v y = -ikp, (11) 
sb x = ikB y v x , (12) 
sb y = -B y d x v x , (13) 
d x b x + ikb y = 0, (14) 

where k is the wavenumber of perturbation and s is the vari- 
able associated with the temporal Laplace transformation. 
Note that these equations describe the perturbations outside 
the forced regions. The effect of the force term is taken into 
account when the perturbations in the different regions are 
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connected across the forced regions (when deriving Eqs. 
and l25l ). Eliminating p, b x , and b v in Eq. ( TTOb by using 
Eqs. dTTb-(fT3T> and (O, we obtain the following ordinary dif- 
ferential equation for v x : 



we get 



(15) 



(4>Dl-xk 2 ) v, = 0, 

(j) = B]k 2 + s 2 (B 2 z +W), 

X = B 2 (k 2 + S 2 ) + S 2 (B 2 + W), 



which has the general solution of the form 

v v = a exp(y / x/ ( l ) kx) + f3 exp(— \Jx/4>kx). (16) 
Other variables are expressed by using v x , e.g., 



p = -s(W+B 2 )d x v x /k 2 , 
b y =-B y d x v x /s. 



(17) 
(18) 



Because the perturbations must vanish at x = ±oo, B x (f3 in 
region 1) and a?, (a in region 1) should be null coefficients. 
We write v x in each region as follows: 



v xA = a x e^ 

v x3 =foe-^ k (*- A \ 



>xl<t>kx 



(19) 
(20) 
(21) 



Let us consider junction conditions for perturbations at the 
interfaces. At the moment, we have four undetermined am- 
plitude of perturbations a\ , a 2 , p\, and j3 3 , indicating that to- 
tal four junction conditions are necessary to determine them. 
Two of them are the continuity of the velocity normal to the 
interfaces that leads to v Xl i(0,i) = v x ^(0,s) and v x ^ 2 (A,s) = 
v Xl 3(A,s) to the first order of the perturbed variables. Sub- 
stituting ([T9ll-(l2T1i to the conditions, we obtain 



ai = a 2 exp(-Vx/#A) + /3 2 , 
a 2 + Pi &Kp(-xfxJ~4>k A) = /3 3 • 



(22) 
(23) 



The remaining two conditions, which describe the force bal- 
ance at each interface, are obtained by integrating (01 with 
regard to z across each inhomogeneous region. After the tem- 
poral Laplace transforming of the integrated equations, to the 
first order of the perturbed variables, we obtain 

(p 2 [0, s]+B ya b ya [0 7 *]) - (p! [0, s] +B yA b yA [0, s]) 

= (w tot , 2 -w tot , 1 )ya,o, (24) 

(p 3 [A, s ]+By 3 byAA,s])-(p 2 [A, s ]-By 2 b > , 2 [A,s]) 

= (W tot , 3 -W tot , 2 )y&.o, (25) 

where £o is the initial amplitude of the interface, and we have 
defined the total enthalpy W to t = W + B 2 +B 2 . In the above 
expression, we have used the fact that v x is continuous across 
the interface and also the case of b x due to (1121 1. Substituting 
(O-dJTJi into (Hill and (|25j via ([T7J» and ([THJ, we obtain two 
conditions written by on,a 2 ,f3 2 , and /J3. Solving these two 
conditions, (l22l . and (l23l with respect to the four coefficients, 



ol\ = ksV(W tot , 2 - 



03 



c 

Aisfi + Bi 



Wtotl), 



A.O 



C 



ksV(W tot . 2 -W tot i), 



(26) 
(27) 



A = {W,ot,2+W toM -2B^+(W to t,2-W toM ) c y s l 
+2B 2 y k 2 . / 

B = {2 (W tot , 2 - B 2 ) s 2 + 2 B 2 k 2 }e~ kA : 

C = {(W tot , 2 + W totil -2B 2 y ) s 2 + 2B) k 2 } 2 

-(W tot . 2 -W totA fe- 2kA s\ 

where we have used the zeroth order conditions W tot ,3 = W tot .i 
and By = B y \ = B y 2 = Byj, and again assumed slow motion of 
the perturbed flows compared to the speed of light: x/4 1 — 1 ■ 
The coefficients a\ and are equivalent to the 
Laplace transform of the interface velocities £[3,^a(0] an d 
£[3,£b(0L respectively. The inverse Laplace transform of 
d26| i and (f27) yield the temporal differential equations for the 
two interfaces: 



at 



x 



A =2{W t 



W tot ,i)kVH(t) 

| B ,o)cos(2 yfjlBykt) 
kA MAfl + ko)cos(2V^B y kt) 



»(l+e- kA )(L 



(28) 



lot. 2 



(1- 



-* A ) + W to u(l+^ A )-2^}, 



-1 



■ 2{ W tot , 2 ( 1 + e* A ) + Wtot, 1 ( 1 " e* A ) - 2 B 2 } , 



where the equation for the interface B is obtained by chang- 
ing the overall sign of the right hand side of ( |28l and by inter- 
changing the subscripts A and B. 

4. PROPERTIES OF SOLUTION 

4.1. Short wave-length limit 

Let us first consider small-scale perturbations compared to 
the thickness of the current sheet, which leads to independent 
evolutionary equations of the two interfaces. Taking the limit 
of k A — > 00 in (l28l and integrating with regard to t , we obtain 

a,B (0 = a,B (0) ± AV kr a,B (0) sin (t/r) , (29) 

Wtot,2-Wtot,l 



A-- 



W toL2 + W totA -2B 2 / 



(30) 



r = Wot,2 + W to t,i-2B?) 



2B y k). (31) 

where plus (minus) sign in the second term of the right hand 
side of ( |29l is for the interface A (B), and r represents the 
lateral Alfven crossing time. Eq. ( 1291 shows that, as long as 
t <^ t (and A ^ 0), interface deformation grows linearly with 
time with the velocity of 



£ = AVk£(0)e 



iky 



(32) 



For t > r, the interface oscillates at the Alfven frequency 

1 /2 

oj ~ B y k/W t0 [ due to the magnetic tension force. In the case 
where the guide field is absent (B y = 0), the interface grows 
perpetually with this speed. In the Newtonian limit (p^S> P 
and B 2 ), A reduces to the Atwood number A = (p 2 -p\)/(p 2 + 
pi), and the growth velocity d32l exactly recover that of the 
original Richtmyer-Meshkov instability (Richtmyer 1960). In 



4 



T. INOUE 



ikV-O.O - -- 

lkV-40 

lkr-S0 

lkV-120 

<W=160 

lkV-200 





FIG. 2. — Evolutionally sequence of the solution j28). Left panel is the 
case of kA = 0.1 and §a/§b = 1. and right panel is the case of kA = 1 and 
Ca/Cb = -1. Other parameters are B y = 0, W tol ,2 = 2, Wt t 1 = 1 and £g = 0.1. 

addition, when B y = 0, the condition A > reduces to the cri- 
terion of Lyubarsky's instability. Thus, the coefficient A can 
be regarded as the generalized Atwood number. It is notewor- 
thy that, differently from the Rayleigh-Taylor type instability 
(including Lyubarsky's instability), the Richtmyer-Meshkov 
type instability can grow even when the Atwood number is 
negative, i.e. the interface B is also unstable. The general- 
ized Atwood number of the current sheet in the relativistically 
magnetized plasma (fil j ^> W\ and P% ^> pi) is A = 1/3, indi- 
cating that the interface is always unstable with respect to the 
shock passage. On the other hand, in the Newtonian limit, the 
current sheet can be unstable provided its density is different 
from the external plasma (in the case of the isothermal Harris 
current sheet, the Atwood number is ^4 = (1 +2f3)~ l , where /3 
is the ratio of the thermal to magnetic pressure in the region 
1). 

4.2. Long wave-length limit 

Next, we consider large-scale perturbations compared to the 
thickness of the current sheet (k A <C 1). The Taylor expan- 
sion of (l28l with regard to A: A results in 



d 2 d i 



Wtot,2-Wtot,l T7 ,2 



(33) 

V^A£A,ocos(f/T)//(f), (34) 



r = (W to t,i-^) 1/2 /(fi^), 

where we have used the fact that the initial amplitudes £a.o 

and £b.o must be equal in this limit. Eq. ( 1331 shows that the 
two interfaces move toward opposite directions. The evolu- 
tionally sequence of the interfaces for the long wavelength 
perturbation is shown in left panel of Fig. [2] The growth in 
opposite phase is not surprising, because the generalized At- 
wood number A at the two interfaces have opposite signs. 

4.3. General properties 

In the general case including kA ~ 1, as the cases in 
both short and long wavelength limits, the following items 
are true: (i) The interfaces can be unstable if the general- 
ized Atwood number A is nonzero, (ii) In t <C r or B y = 0, 
where t ~ y/W tot /(B y k) is the lateral Alfven crossing time, 
the interfaces grow with constant speed of ~ AkV£o. (iii) 
In t > t, the interfaces oscillate with the frequency ~ 1/r. 



Interestingly, eq. ( 1281 ) suggests that, for B v = (or t <C t), 
the growth speed of the interface A {d^/dt) becomes zero 

for the special case of £a,o/£b,o = {(W , to t,2 + Wiot,i) + (Wtot ) 2- 
W tot ^)e- 2kA }/(2W tot 2e- kA ) = r c and that of the interface B 

becomes zero for £a,o/£b.o = '"c > wnere r c is always larger 
than unity. Furthermore, when the ratio of the initial am- 
plitudes is in the range r c > £a.o/£b,o > r~ l , the ratio of the 
growth speeds of the interface A and B is negative (out of 
phase growth). When the initial ratio is larger than r c or 
smaller than r~ l including negative value, the ratio of the 
growth speeds is positive (in phase growth). The reason for 
this as follows: Since the generalized Atwood number A has 
opposite signs at the two interfaces, the interfaces at the same 
y are accelerated toward the opposite (same) directions when 
the initial interfaces are in phase (out of phase). However, 
when the magnitude of the initial amplitude is very differ- 
ent, the motion of the interface with minor initial amplitude is 
dragged toward major one, because the volume of the accel- 
erated regions depends on the initial amplitude (see Fig. [T]). 
The evolutionally sequences of the in phase and the out of 
phase growths are shown in the left and right panels of Fig 
respectively. 

5. DISCUSSION 

Finally, we discuss implications of the instability. Re- 
cent particle-in-cell simulation of a relativistically magne- 
tized pulsar wind have shown that shock propagation through 
the stripes of opposite magnetic field polarity induces driven 
magnetic reconnection (Sironi & Spitokovsky 2011). The 
triggering mechanism of the magnetic reconnection could be 
the instability studied in this article. 

In addition, recent MHD simulations have shown that the 
growth of Richtmyer-Meshkov instability induced by a (rel- 
ativistic) shock propagation through inhomogeneous den- 
sity medium generates turbulence in its nonlinear stage (Gi- 
acalone & Jokipii 2007, Inoue et al. 2009, 201 1, 2012, Beres- 
nyak et al. 2009). Because the instability of the relativistic 
current sheet found in this study is also a Richtmyer-Meshkov 
type instability, it is quite reasonable to expect the develop- 
ment of turbulence. 

The instability evolves into the nonlinear stage, if the ampli- 
tude of the interface can grow to £ ~ l/k before the magnetic 
tension force begins to suppress the growth (t < r). Using the 
growth velocity d32l . the above condition is reduced to 



V 

VZy 



> 



(AkHuT 



(35) 



where Va, v = B y j -\AVtot is the lateral Alfven velocity. Note 
that if there is no guide field (B y = 0), which is plausible at 
least in the pulsar wind, the instability can always go into its 
nonlinear stage. As in the case of the Richtmyer-Meshkov 
instability, the growth of the instability is not exponential with 
time but linearly with time (Richtmyer 1960), indicating that 
the timescale of the evolution depends on the initial amplitude 
of perturbation. If we consider the initial perturbation of scale 
k > A -1 with amplitude £b,o = £a,o = £o in which the instability 
is the most influential, the instability grows to nonlinear (£ k ~ 
1) after the following elapsed time from the shock passage: 



(AV&k 1 ) 



2\-l 



(36) 



where we have estimated the timescale based on eq (|28l . 
This suggests that, if the initial amplitude is much smaller 
than the wavelength of perturbation (A = k~ x <C £o)> it takes 
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huge time for the instability to evolve into the nonlinear stage. 
Fortunately, however, we can expect large initial amplitude 
of perturbation owing to the tearing mode instability and/or 
Lyubarsky's instability that can grow prior to the shock pas- 
sage. It is discussed in Lyubarski (201 1) that Lyubarsky's in- 
stability can grow in pulsar winds, gamma-ray burst jets, and 
active galactic nuclei at least in the scales A < A, although it 
may not fully dissipate the magnetic field especially in pulsar 
winds. Thus, we can expect excitation of turbulence quickly 
(on the order of crossing time A V" 1 for £o ~ A) after the shock 
passage. Note that, in Lyubarsky's instability, either the in- 
terface A or B can be unstable depending on the orientation 
of acceleration in the most unstable scales of A < A, while 
our instability can grow in both interfaces. Therefore, even 
if Lyubarsky's instability is going on, our instability that agi- 
tates the current sheet as a hole is necessary in order to excite 
"turbulence". 

Once the current sheet becomes turbulent, the induction 
of so called "turbulent reconnection" (Lazarian & Vishniac 
1999, Kowal et al. 2009) and the effect of turbulent stretching 
of the current sheet (Takamoto et al. 2012) dissipate the mag- 



netic field rapidly regardless of plasma conductivity. The ex- 
cited turbulent reconnection would then evolve through gen- 
eration of more turbulence by the positive feedback of recon- 
nection outflows. Therefore, the instability found in this study 
can be applied as a triggering mechanism of rapid magnetic 
energy release in variety of high-energy astrophysical phe- 
nomena such as pulsar wind nebulae, gamma-ray bursts, and 
active galactic nuclei, where the shock wave is supposed to 
play a crucial role. Differently from the secular instabilities 
such as the tearing mode instability and Lyubarsky's insta- 
bility, sudden onset of our instability due to the shock pas- 
sage and following turbulent reconnection may be preferred 
especially in the intermittent phenomena such as gamma-ray 
bursts. 
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